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Abstract. FEARLESS (Fluid mEchanics with Adaptively Refined Large Eddy Simulations) is a new numerical scheme 
arising from the combined use of subgrid scale (SGS) model for turbulence at the unresolved length scales and adaptive 
mesh refinement (AMR) for resolving the large scales. This tool is especially suitable for the study of turbulent flows in 
strongly clumped media. In this contribution, the main features of FEARLESS are briefly outlined. We then summarize the 
main results of FEARLESS cosmological simulations of galaxy cluster evolution. In clusters, the production of turbulence is 
closely correlated with merger events; for minor mergers, we find that turbulent dissipation affects the cluster energy budget 
only locally. The level of entropy in the cluster core is enhanced in FEARLESS simulations, in accord with a better modeling 
of the unresolved flow, and with its feedback on the resolved mixing in the ICM. 
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INTRODUCTION 

The formation and evolution of the large scale structure, which proceeds through the process of hierarchical clustering, 
plays a key role for the conversion of potential gravitational energy into kinetic and internal energy in galaxy clusters 
and groups. In this process, cluster mergers induce bulk motions in the intra-cluster medium (henceforth ICM), with 
velocities of the order of 1000 km s~ 1 . The shearing instabilities associated with the merger events and the ubiquitous 
shock waves in the ICM are thus expected to stir the gas and make the flow turbulent. 

Cluster formation is a typical example of turbulence generation in a strongly clumped medium. As the Eulerian 
approach is generally considered to be superior to the Lagrangian one in modeling hydrodynamical instabilities 
HI I2i 1211), grid-based numerical codes are a suitable tool to address this problem, especially when the Adaptive Mesh 
Refinement (henceforth AMR, 0>H,|^]) technique is used. On the other hand, as in most of the astrophysical fluids, 
the Reynolds number in the ICM is estimated to be Re 1 (although with many issues on the value of the kinematic 
viscosity; 0, SEE Ell EH] X implying a number of degrees of freedom of the dynamical system lfl2j[l3ll which is not 
manageable even in state-of-the-art numerical simulation^ Therefore, even with AMR, it is usually not possible to 
resolve all the turbulent cascade, down to the dissipative length scale foil . 

In many fields of computational fluid dynamics, the influence of unresolved turbulence on the resolved scales is 
modeled by means of heuristic subgrid-scale (SGS) models, coupled to the large scales of the system, for which the 
fluid equations are solved (Large Eddy Simulations, LES; 11410 . LES are customary, for example, in simulations of 
Type la supernovae lfL^[l^[T^[l8l[l9tl . Recently, ||20il developed a similar tool designed for studying Rayleigh-Taylor 
driven turbulence in the interaction between AGN outflows and the ICM lf2lll22tl . 

In this contribution, we describe a novel numerical tool that unites LES and AMR, called FEARLESS (Fluid 
mEchanics with Adaptively Refined Large Eddy Simulations). FEARLESS combines the adaptive refinement of the 
regions where turbulent flows develop with a consistent modeling of the SGS turbulence, and is argued to be a suitable 
tool in simulations of turbulent clumped flows. The numerical formalism and its first application to the physics of 
galaxy clusters is presented in l23ll . The main features and results are collected here. 



In case of cosmological simulations, the problem is even more severe because the evolution of galaxy clusters has to be followed in the 
cosmological context of their collapse, i.e. in simulated computational domains much larger than the cluster scale itself. 



NUMERICAL TOOLS 



FEARLESS has been implemented into the ENZO cod^E v. 1 .0 |l24l] . an AMR, grid-based hybrid (hydrodynamics plus 
N-Body) code based on the PPM solver l25ll modified for the study of cosmology [26] , which provides the necessary 
infrastructure for performing numerical simulations of cosmological structure formation. 



Adaptive mesh refinement of turbulent flows 

As a first step in the development of FEARLESS, some work has been devoted to the development of refinement 
criteria which are best suited for refining turbulent flows. One such criterion was introduced in M27I1 . and is based on 
the regional variability of structural invariants of the flow, i.e. variables related to the spatial derivatives of the flow 
velocity. An example is the modulus of the vorticity 0) = V x v (the curl of the velocity field), expected to become 
high in regions where the flow is turbulent. The regional threshold for triggering the refinement is expressed in the 
comparison of the cell value of the variable q(x,t) and the average and the standard deviation of q, calculated on a 
local grid patch: 

q{x,t)>{q)i{t) + aXi{t) (1) 

where A, is the maximum between the average (q) and the standard deviation of q in the grid patch i, and a is a tunable 
parameter. 

This technique has been profitably used in simulations of idealized minor mergers i28ll and in full cosmological 
simulations of galaxy clusters, focused on minor merger events 12911 . In both cases, the new AMR criteria provide a 
better resolution of the turbulent flow, as demonstrated by larger values of the velocity dispersion on resolved scales 
and (for l29ll ) by an increased floor of entropy in the cluster core, due to an increased mixing. 

The new criteria have been used also in a recent study of major mergers 13QL allowing to follow the evolution of 
merger shocks and the related post-shock region. The evolution of turbulence in that region and in the cluster core has 
been explored, and interesting morphological comparisons with the observed symmetrical radio relics in the merging 
cluster Abell 3376 l3~lll have been performed. 



Subgrid scale model 

Here, we give a general and qualitative account of the main features of the SGS model; we refer to lfl3ll and l23ll for 
a more detailed and rigorous description. 

As shown by lfl3ll . the equations of motion for a compressible, viscous, self-gravitating fluid can be decomposed 
into large-scale (resolved) and small-scale (unresolved) parts, using the Germano filter formalism ll33ll . By means of 
filtering, any field quantity a can be split into a smoothed part (a) and a fluctuating part a', where (a) varies only at 
scales greater than the prescribed filter length. 

Following this procedure, one can define the SGS turbulent energy e t , which can be interpreted as a subgrid energy 
buffer between the resolved kinetic energy and the internal energy (graphically shown in Fig.Q]). The turbulent energy 
is governed by an equation of the following form: 

^- t {p)e t + ^-vj(p)e t = D + E + r-(p)(A + e), (2) 

where p is the gas density and vj is the density- weighted filtered velocity component, according to Favre |34{j. The 
quantities on the right-hand side of Eq. [2] determine the evolution of e t and are the turbulent diffusion term D, the 
turbulent production term E, the pressure dilatation term A and the viscous dissipation term e. Their role in the energy 
budget is visually described in Fig.Q] and the way they are modelled (i.e. their so-called closures) represent the core of 
the SGS model. The filtering procedure couples the hydrodynamical equations to the unresolved scales by the inclusion 
of terms from Eq.|2] as described in 1IT31 12311 . 



ENZO homepage: http://lca.ucsd.edu /so ft ware /enzo/ 
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FIGURE 1. Graphical description of the energy budget in FEARLESS, showing the energy components and the exchange terms 
between them. From 13211. 



Combining AMR and LES 

A severe limitation of the SGS model is the implicit assumption of the isotropy of the velocity fluctuations on 
subgrid scales. This limits the LES methodology to flows where all anisotropics stemming from large scale features, 
like boundary conditions or external forces, can be resolved. In FEARLESS, this issue is cured by the adaptive 
local adjustment of the grid resolution, in order to ensure that the anisotropic, energy-containing scales are resolved 
everywhere as much as possible. In this way, it is assumed that turbulence is asymptotically isotropic on length scales 
comparable to or less than the grid resolution. 

When a grid location is flagged for refinement in Enzo, a new finer grid is created, and the cell values on the finer 
grid are generated by interpolating them from the coarser grid using a conservative interpolation scheme. At each 
timestep of the coarse grid, the values from the fine grid are averaged and the values computed on the coarse grid 
(in the region where fine and coarse grid overlap) are replaced. An additional step of the refinement procedure in 
FEARLESS is a scale correction in the energy budget, resulting in a energy transfer from the SGS turbulent to the 
resolved kinetic energy (or the opposite for grid derefinement). The relation between the SGS turbulent energies at 
different levels of refinement l\ and 1% is given by the assumption of Kolmogorov scaling l35ll36ll : 



et.2 \h) 
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(3) 



The required correction steps for the resolved velocity components and for the turbulent velocity in the newly 
created fine grid (the primed quantities in the following equations) are: 



e t = e t r A (5) 

where e^ n is the resolved kinetic energy and r\ is the refinement factor of the mesh. The resolved energy is adjusted 
such that the sum of resolved energy and turbulent energy remains conserved. The opposite procedure is applied at 
grid derefinement. 



GALAXY CLUSTER SIMULATIONS WITH FEARLESS 



The numerical technique presented in the last section has been applied to cosmological simulations of the formation 
and evolution of a galaxy cluster ll23ll . Two cluster simulations, run with and without the use of the SGS model but 
otherwise with an identical numerical setup, are compared. 

A fiat ACDM background cosmology is assumed, with parameters £2a = 0.7, £2 m = 0.3, D.b = 0.04, h = 0.7, 
Og = 0.9, and n = 1. The simulations have been initialized at redshift Zm = 60 using the transfer function by [37], 
and then evolved to z = 0. 

The computational box has a size of 128 Mpc h-\ The root grid (level / = 0) has 128 3 cells and 128 3 N-body 
particles. A static child grid (/ = 1) is nested inside the root grid with a size of 64 Mpc h , 128 3 cells and 128 3 N- 
body particles. The mass of each particle in this grid is 9 x 1O 9 M0 hT 1 . Inside this grid, AMR is allowed in a volume 
of 38. 4 Mpc ft -1 , using the overdensity refinement criterion as described in l29ll with an overdensity factor / = 4.0. 
The refinement factor between two levels was set to 2, allowing for an effective resolution of 7.8 kpc h~ l . 

The static and dynamically refined grids were nested at the place of formation of a galaxy cluster, previously 
identified in a low-resolution run using the HOP algorithm ll3~8ll . The cluster did not experience in its recent history 
(z < 1) any major merger, therefore it was considered an ideal case to study the role of minor mergers for injecting 
turbulent energy in the ICM. 

In Fig. [2] a short time series of density and SGS turbulent velocity slices is presented. Several minor mergers in the 
ICM can be identified, especially at the locations where the SGS turbulence is large. Indeed, a considerable amount 
of turbulent energy is localized in front and in the wake of the merging clumps (panels[2j? and d). The morphological 
evolution in Fig. [2] provides therefore a remarkable insight on the markedly local behavior of the production and 
dissipation of turbulence in clusters. 

The energy content of the turbulent flow in the ICM is rather small, in agreement with the subsonic nature of the 
flow. The SGS turbulent energy is generally smaller than 1 % of the internal energy, while the energy associated with 
the velocity dispersion at the length scale of 100 kpc h is at the level of a few percent. On the other hand, the 
turbulent contribution to the energy budget is more relevant at the location and in the wake of the merging subclumps. 
According to l23ll . the value of the SGS turbulent energy in the vicinity of a subclump is one order of magnitude larger 
than the typical average value in the ICM. In addition to the numerical dissipation of the code, which converts kinetic 
energy directly into internal energy, there is an additional dissipation channel due to the energy flux from resolved 
scales to the internal energy through the SGS energy buffer. 

An interesting result, closely linked with the point just discussed above, comes from the comparison between radial 
profiles in the FEARLESS simulation and in the standard adiabatic run (Fig. [3). A change in the temperature and 
density profiles at the cluster center is clearly visible. In particular, in the FEARLESS run T is larger (3%) for central 
distances r < 0.07 R v i r , with respect to the standard run (Fig. 0?). Consequently, the gas in the core is less dense 
(Fig. [3J;), so that the ICM locally remains in hydrostatic equilibrium. The local energy budget in the cluster core is 
therefore modified by the SGS model. 

The effects of the SGS model is reflected in the entropy which is defined, as customary in cluster physics, as 



where 7=5/3 and p is the gas density. The entropy in the cluster core is higher in the FEARLESS run as compared 
to the standard run (Fig. |3J:). This result is produced by the combination of two factors: the enhanced local dissipation 
of SGS turbulent to internal energy, provided by the SGS model, and a moderate feedback of the SGS model on the 
resolved scales, increasing the magnitude of the resolved mixing (cf. I23I1 for more details). 

The enhanced entropy level in the core is in agreement with known results of the comparisons between SPH and 
grid codes on this issue l39l l40l l4lll . Whenever the modeling of the turbulent flow in cluster simulations is improved, 
either by using grid codes rather than SPH, or by especially designed AMR criteria for refining turbulence 12911 . or by 
using a SGS model i23"ll . the level of core entropy is increased. Recently, l42ll pointed out that the core temperature 
and entropy in grid-based codes are affected by a spurious increase, caused by the N-body noise in the gravitational 
force field. In our opinion, further investigations are certainly needed to disentangle physical from numerical effects 
on this issue. 




FIGURE 2. Slices of baryon density (left-hand panels, a and c) and turbulent velocity q = yf2e{ (right-hand panels, b and d) at 
different redshifts z, for the cosmological simulation run with FEARLESS. The density is logarithmically color coded, whereas q 
is linearly coded in km s , according to the colorbar on the left of panel b. The overlayed contours show density. The slices show 
a region of 6.4 x 6.4Mpc h around the center of the main cluster followed in the simulation. Panels a and b refer to z = 0.05, 
panels c and d to z = 0. 

SUMMARY AND CONCLUSIONS 

Large Eddy Simulations (LES) are based on the notion of filtering the fluid dynamic equations at a specific length 
scale. A separation between the resolved and the unresolved length scales of the flow is thus performed. The latter 
are treated by means of a subgrid scale model, which is coupled to the hydrodynamical equations governing the 
former. In principle, a single scale separation is not practical in simulations of clumped media, and moreover it is 
not immediately compatible with the concept of adaptive mesh refinement (AMR), often used to study astrophysical 
systems. The solution proposed to this problem is the development of a new numerical scheme exploiting AMR and 
LES in combination, called FEARLESS. 

In this contribution we briefly reviewed the basic ideas of the FEARLESS approach (AMR, SGS model and their 
coupling), and then the application of this technique to galaxy cluster simulations was discussed. The subsonic flow in 
the ICM are do not have a dominant role in the energy budget (as confirmed by previous analyses, e.g. J43il39ll44ll45ll ). 
and therefore the impact of the SGS energy content is globally negligible. Nevertheless, at the location of the merging 
subclumps, the role of turbulent dissipation become sizeable. Such local effects are particularly relevant in the cluster 
core, where an enhanced entropy floor results in the FEARLESS simulation. 




FIGURE 3. Cluster radial profiles of selected quantities at z = 0. The dotted line refers to the simulation without SGS model, 
whereas the solid line is for the simulation using FEARLESS. The radii on the x-axis are expressed in fraction of the virial radius, 
R v i r = 1.37 Mpc h~ l for the considered cluster. Panel a reports the gas density, in b is the mass-weighted temperature and in c the 
gas entropy, defined in the text. 



Future work in this area will regard cosmological simulations of major mergers. These events are particularly 
interesting because the turbulent energy injected by them has a significant role in the energy budget, with respect 
to minor mergers. As a first step, a study which only makes use of AMR, triggered by the regional variability of the 
compression rate of the flow, has been performed 113011 . Furthermore, it will be also important to explore the turbulent 
flow not only in the ICM, but also in the cluster outskirts and in the cosmic web. The knowledge of the turbulent 
state of this gas is potentially important for shaping the emission and absorption lines associated with the Warm- 
Hot Intergalactic Medium (see 14611 f or a review). From a computational viewpoint, numerical improvements in the 
FEARLESS scheme are ongoing 14711 . in order to better account for complex flows with large density gradients. 
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